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Abstract 

An important ingredient for the non-equilibrium evolution of partons at RHIC and LHC is to 
have some physically reasonable initial conditions for the single particle phase space distribution 
functions for the partons. We consider several plausible parametrizations of initial conditions for 
the single particle distribution function fi(x,p) and fix the parameters by matching / f(x,p)p fl da fl 
to the invariant momentum space semi-hard parton distributions obtained using perturbative QCD 
(pQCD), as well as fitting low order moments of the distribution function. We consider parametriza- 
tions of fi(x,p) with both boost invariant and boost non-invariant assumptions. We determine the 
initial number density energy density and the corresponding (effective) temperature of the minijet 
plasma at RHIC and LHC energies. For a boost non-invariant minijet phase-space distribution 
function we obtain ~ 30(140) /fm 3 as the initial number density ~ 50(520) GeV/fm 3 as the 
initial energy density and ~ 520(930) MeV as the corresponding initial effective temperature at 
RHIC (LHC). 

PACS numbers: PACS: 12.38.Mh; 12.38.Bx; 13.87.Ce; 25.75.-q 
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In order to study thermalization of quark-gluon plasma at RHIC and LHC in any kinetic 
approach, one needs suitable initial conditions for the single particle phase-space distribution 
of partons. The hard and semi hard parton production can be calculated by using perturba- 
tive QCD (pQCD) [lj]. The soft parton production cannot be calculated in a kinetic approach 
within perturbative QCD. One approximate strategy for looking at soft parton production is 
to assume the creation of a classical chromofield 0,0,0]. The first problem one faces in spec- 
ifying initial conditions for the phase space single particle distribution function f(x,p,to) of 
partons is that from pQCD one has information only about the momentum space distribu- 
tion of minijets. The pQCD calculation (using scattering matrix element squared and the 
structure functions) is done in an approach which calculates probabilities between initial 
and final states so that any information about space-time evolution is lost. Thus we need 
to supplement pQCD with assumptions about the dependence of f(x,p) on x. 

In this paper our strategy will be to use pQCD to determine the single particle distribution 
function for the partons in momentum space, and then by taking simple parametric forms for 
f(x,p) constrain the parametrization by fitting the single particle distribution as well as the 
average transverse momentum. This process can be refined by adding more parameters and 
fitting further moments of the jet distribution function. Of course, one hopes that the details 
of the parametrization of the initial conditions do not significantly effect the thermalization 
process. 

The lowest order pQCD inclusive (2 — ► 2) minijet cross section per nucleon in A-A 
collision is given by: 



where &ij~>ki is the elementary pQCD parton cross sections for the process ij — > hi. As 
gluons are the dominant part of the total minijet production, we will only consider the 
processes gq(q) —>■ gq{q) and gg — > gg in this paper. The partonic level cross sections for 
these processes are given by: 
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with the kinematic relations: 

s = X1X2S, and i — — -[1 — tanh( ^ 1 - — )]. (5) 

To compute the minijet cross section using Eq. we need to specify the minimum trans- 
verse momentum cut-off po above which the incoherent parton picture is applicable. The 
momentum cut-off is ~ 1 GeV at RHIC and ~ 2 GeV at LHC, which are obtained from 
the argum e„t that the, b Nation of the gl „o„ structure fuuetiou at ,w , Q For 
our quantitative calculation we use po= 1 GeV at RHIC and 2 GeV at LHC. The nuclear 
modified parton distribution functions which includes the shadowing effects are given by 
fi/A(x,Q 2 ) = RA(x,Q 2 )fi/N(x,Q 2 ) where A and N stand for nucleus and free nucleon re- 
spectively. In this paper we use the GRV98 set of parton distributions for the free nucleon 



distribution function fi/N(x,Q 2 ) f| (in momentum space). The GRV98 analysis uses the 
low x HERA data on deep inelatstic scatterings along with data from other hard scattering 
processes at fixed Q. The Q 2 evolution of the parton distribution function is performed by 
using perturbative QCD evolution equations. We use the EKS98 numerical parametrizations 
for the ratio function Ra{x, Q 2 ) [7[. The EKS98 parametrization uses NMC and E665 struc- 
ture functions from deep inelastic lepton-nucleus collisions and E772 Drell-Yan data from 
proton-nucleus collisions at fixed Q. The Q 2 evolution of the structure function is studied 
by using the DGLAP evolution equation. We multiply the above minijet cross section by a 
standard K factor (K=2), to account for the higher order contributions. 

For central collisions, the minijet cross section (Eq. (£[])) can be related to the total number 
of partons (N^ et ) by 

dNi et T ^ f 2np T 

ijkl 

where T(0) = 9 A 2 /8irR A fm is the nuclear geometrical factor for head-on AA collisions (for 
a nucleus with a sharp surface). Here Ra = 1.L4 1 ' 3 is the nuclear radius. Similarly the 
transverse energy distribution of the minijets are given by: 

dE J T et r ^. n . r 2np 2 



K T(0) J rf|/2^Ell fi/A{Xl,PT) x 2 fj/A{x2,PT) °ij->kl{s , t, u) , (6) 



dydpi 



/ 27rp 
dy2—^J2 x i fi/A{x\,p 2 T ) x 2 fj/A{x2,P 2 T ) Vij->ki(s,i,u). (7) 

S ijkl 

Our numerical results for the pt distribution of the minijets computed from eq. © are 
plotted in Fig. 1. 
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FIG. 1: Transverse momentum distribution of the initial minijets at RHIC and LHC 

The rapidity distributions of the minijet production at RHIC and LHC are obtained 
by integrating over pt for px > Po are shown in Fig. 2. It can be seen that the minijet 
distribution is not flat over the whole rapidity range. Hence it is not a particularly good 
approximation to assume boost invariance |s| and restrict oneself to the mid rapidity region. 
One important question that the transport equations should answer is how different the final 
rapidity distribution of the hadrons (or other signatures) from the rapidity distribution of 
the minijets at the pre-equilibrium stage are Of course this needs to be supplemented 
by further dynamical evolution near the hadronic phase transition. For this reason it is 
not appropriate to take into account just the flat minijet rapidity distribution seen in the 
mid rapidity region when evolving the minijet plasma. Partons formed outside the mid 
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FIG. 2: Rapidity distribution of the initial minijets at RHIC and LHC 

rapidity region do propagate in the pre-equilibrium stage with their rescatterings and hence 
one must consider all the partons in determining the initial condition of the QGP in the 
pre-equilibrium stage. 

Thus we require additional information about how these partons are distributed in coor- 
dinate space. A simple ansatz based on scaling ( 8]) allows extracting energy density from 
the transverse energy Et distribution via: 

ttR z a t dy 

This formula is correct if partons are uniformly distributed in coordinate space and the mo- 
mentum rapidity y is equal to the coordinate rapidity 77. For an 1+1 dimensional expanding 
system coordinate rapidity 77 is defined via: t = r cosh 77 and z = rsinhr/. 



Traditionally the above formula is widely applied to obtain energy density from Et dis- 
tribution of minijets. Although the above formula might be applicable in a situation where 
the minijet distribution is flat in rapidity, it is precisely the scale breaking effects that we 
hope to understand from our transport approach. As it can be seen from Fig. 2 the initial 
minijet distribution function is not flat in the entire rapidity range and hence the above 
formula is not applicable in the pre-equilibrium region. Also it is not clear that partons 
are uniformly distributed at the initial time and just dividing by a volume (see eq. (|5jl) to 
obtain the energy density from the total energy is a rather crude estimate of initial condi- 
tions. To get a better estimate for the initial conditions one should take into account the 
correlation between coordinate rapidity 77 and momentum rapidity y. We will assume this 
correlation has a particular form, and then look at the sensitivity of the global quantities to 
the parametrization. 

For an expanding system the mini jet number distribution eq. © can be related to the 
phase-space distribution function via 
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d 3 N jet f 

^a? T =scjf^ P )p"d^ ( 9) 

where gc = 16 is the product of spin and color degrees of freedom, 

da^ = n R 2 A rdr] (cosh r], 0, 0, sinhr?), (10) 

and 

P^ — (pt cosh y,pr cos 4>,Pt sin (p,pTsmhy), (11) 

for an 1+1 dimension expanding system. Using the above relations we get at the initial time 
To (=l/po): 

dN jet r 

j j 2 = 9cnR 2 A T / dr] p T cosh(r7 - y) f(p T , 77, y, r ), (12) 
ndydp T J 

where we assume a transverse isotropy at the early stage. In the following we will obtain the 
initial phase-space distribution function of the gluon minijets by using both boost invariant 
and boost non-invariant schemes. 



A. Boost Invariant Initial Distribution Function From Minijets 



First let us assume that a boost invariant description [8( is appropriate. Then the gluon 
distribution function depends only on the boost invariant parameters r, p T and £ = rj — y. 



We will parametrize the r] — y correlation as Gaussian: 
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f{PT,V-y,r ) = f{p T )e (13) 

where f(pr) wm be determined from eq. (|12|) by using eq. (JHJ) which is obtained by using 
pQCD. Using the above equation in eq. (|T2"j) we get: 

dN^ et r r°° _(v-yf 

= g c 2ir 2 R 2 A T p 2 T f(p T ) J dy J ^ drj cosh{r) - y)e (14) 

which gives: 

— = g c 2n 2 ^R 2 A T p 2 T f(p T )a(p T )e^)/* J d y (15) 

where the maximum allowed value of the momentum rapidity at RHIC and LHC depends 
on px and \fs. Integrating over y we get from the above equation: 

f(Pr) = ^ —, - s , (16) 

g c 2^^R 2 A T p 2 T 2ln(^/2p T + Js/Ap 2 T - 1) *(Pr) 



which gives a boost invariant initial distribution function: 



g <T 2 (p T ) 



/(PT ' ^ g c 2^^R 2 A r,p\2 ln(^/2p T + Js/4p 2 T - 1) ^e^/ 4 " 



Using this boost invariant initial phase-space distribution function of the gluon the initial 
minijet number density is given by: 

n(r ) = 9c J dTp^u^fipTiV - Vi T o) = 9c J d 2 p T p T J G?£cosh£,f(p T , £, r) 

dNi et 

= 27T / dp T (18) 

J 2tc 2 R 2 a t 2 \n( v r s / 2 p T + Js/4p 2 T - 1) 



where dT = ^f. Similarly the initial energy density is given by: 

e(r ) = 9c J dT(p IJ 'u tl ) 2 f(pT,V ~ V^o) = 9c J d 2 p T p 2 T J d£ cosh 2 £ f{p T , £, r) 



dNi e 



2tt / dp T p T ^ ( e -* 2 (pr)/4 + e 3ff2( ^ )/4 ). (19) 

^ 47T 2 i%T 2 ln(v^/2p r + - 1) 



B. Boost Non-Invariant Initial Distribution Function From Minijets 

In the above we used a boost invariant distribution function f(r, £,Pt) as the initial 
phase-space distribution function. However this is not consistent with the actual rapidity 



7 



distribution function obtained from pQCD. To relax the boost invariant assumption, we 
will introduce a function f(pr,y) instead of /(pr) in eq. (|T5jl. The form of the boost 
non-invariant distribution we take is: 

fiPr, V, V, to) = fipr, y)e °^t) , (20) 

where f(pr,y) is obtained from ■ Note that f(pr) in eq. (fT3|l is obtained from 
after integrating over y. However, in the above form (eq. ()2()jl ). the y dependence is more 
general which for small a 2 is close to that obtained by using pQCD. Using the above form 
in eq. (fT2"j) we get: 

dW et 
dydp T 

which gives: 



/oo _ {v~y> 

dr)cosh(ri-y)e ^r) (21) 
-oo 



(/ A " ' gc27r 2 ^R 2 A r p 2 T f(p T , y)a{p T )e^^\ (22) 



dydpT 

From the above equation we get: 



/(Pr ' V) g c 2n^R 2 JopXcr(p T )e° 2 M/± (23) 
which gives a boost non-invariant initial phase-space gluon distribution function: 

Using the above boost non-invariant gluon distribution function the initial minijet number 
density is given by: 

n(r , v) = 9c J dTp^u^fipr, rj, y, r ) = 9c J d 2 p T p T J dycosh(rj - y)f(p T , V, y, T o) 

dw<* 



= 2n j "dp T j 'dy 27r2 ^ irQ .J^^L)^ cosh ^ - Vl (25) 
and the initial minijet energy density given by: 

e( T o,v) = 9c J rfr(p%) 2 /(PT,^,l/,T-o) =9c j d 2 p T p 2 T J dy cosh 2 (r/ - y)f(p T , V, y, T o) 



dW- — 2 



(rj-VT 

271 1 dpTPT I dy 2,\m 2 r n a(i)X^ C ° Sh2{v ~ V) - (26) 



2n 2 ^R 2 A T a(j> T )e° 2 ^)/± 
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C. Initial Energy Density and Number Density at RHIC and LHC 



Using the above expressions for the minijet initial phase-space distribution functions we 
can now predict the initial energy density and number density of the minijet plasma at RHIC 
and LHC As temperature can not be defined in a non-equilibrium situation we define an 
effective temperature by: 

T eff = [^] 1/4 ; (27) 
with the non-equilibrium energy density obtained above. 

Before we calculate the initial conditions we would like to find out the values of the 
unknown parameter a 2 appearing in the above equation. We will fix this unknown parameter 
by equating the first pt moment of the distribution function with the pQCD predicted Et 
distribution as given by eq. (|7Jl. We have for the boost invariant case: 



/ 



JF jet r r Vt — 

r _0 2 D 2 I j 3 f( \ I J F1 dp T 



dp T ^- = 2g C 7r 2 R 2 A T / dp T p s T f(p T ) = / dp T . apT & 

dPT J J v 0F21n(v^/2p T + Js/Ap 2 T - l)a{p T )e° 2 ^)/¥ 



and for the boost non-invariant case: 

dE jet Pt dNJet 

I dPT I dy djv^ = 2 ^ 2 R\t, J dp T J dy p* T f(p T , y)=J dp T J d y -^Z^_ (2 9) 

We determine the value of a by using the pQCD values of , J and , T from eq. 
(JZj) in the left hand side of the above equations. Assuming a 2 independent of px we 
find from the above equations: a 2 = 0.0034(0.0015) at RHIC(LHC) for boost invari- 
ant case and a 2 = 0.28(0.28) at RHIC(LHC) for boost non-invariant case. Note that 
for boost non-invariant case a 2 is same at RHIC and LHC as the normalization factor 
y / 7rcx(pr)e cr2< - PT ^ //4 in eq. (}29|) is energy independent whereas the corresponding normalization 
factor ypUl h\{s/s/2pT + J s/Ap\ — l)a(pT)e a2( - PT ^ 4: in eq. depends on energy. Using the 



above values for the boost invariant case in eq. we find eo=31 GeV/fm 3 with T e ff— 
460 MeV at RHIC and e =287 GeV/fm 3 with T eff = 800 MeV at LHC. The number density 
for the boost invariant case obtained from eq. (JTSJl is 72 =20/fm 3 at RHIC and 85/fm 3 at 
LHC. Using the values of the a 2 for boost non-invariant case we find from eq. eo=51 
GeV/fm 3 with T e// =521 MeV at RHIC and e =517 GeV/fm 3 with T e// =932 MeV at LHC 
for rj=0. The number density for the boost non-invariant case obtained from eq. (|23|) is 
no=30/fm 3 at RHIC and 138/fm 3 at LHC for r]=0. It can be noted that the energy density 



9 



for the boost invariant case is less than that of the boost non-invariant case. This is because 
we have assumed a boost invariance in the entire rapidity range. However, the actual pQCD 
minijet rapidity distribution is not flat in the entire rapidity range as can be seen from Fig. 
2. Therefore, in order to maintain the boost invariance in the entire rapidity range the den- 
sity has to be smaller. However, in boost non-invariant case the rapidity distribution of the 
minijet is minimaly altered and the pQCD form of d ^ dy is used throughout the calculation. 

Let us mention briefly other approaches for obtaining initial condition at RHIC and LHC. 
Commonly used way is to obtain ene rgy density from eq. (JBJ) by using pQCD estimate of 
the minijet Et distribution ^, 11, 12, lj|. This formula assumes that the initial partons are 
uniformly distributed in the coordinate space so that one divides the volume to obtain energy 
density. In [14] a boost invariant Boltzmann form f(x,p,t ) = e~ PT cos H r i-y)/ T was used for 
the initial distribution function with rj — y (boost invariant quantity) correlation taken into 
account. In the parton cascade model 15[ momentum space structure function along with 
coordinate distribution of the nucleons inside the nucleus was used directly in the transport 
equation which is different from the minijet calculation by using pQCD. In the approach 
of this paper we have obtained a phase-space initial non-equilibrium distribution function 
f(x,p) for the boost non-invariant case with minimal altering the and y distribution of 
the minijet partons obtained by using pQCD calculations. 

To conclude, we have obtained a physically reasonable initial phase-space distribution 
function of the partons formed at the very early stage of the heavy-ion collisions at RHIC 
and LHC which can be used to study equilibration of the quark-gluon plasma 
llfi| . To obtain such phase-space distribution function we have considered several plausible 
parametrizations of fi{x,p) and have fixed the parameters by matching / f(x,p)p fl da f j i to 
the invariant momentum space semi-hard parton distributions obtained using QCD (pQCD) 
as well as fitting low order moments of the distribution. Once the parameters of f(x,p) 
are found, we then have determined the initial number density, energy density and the 
corresponding (effective) temperature of the minijet plasma at RHIC and LHC energies. 
For a boost non- invariant minijet phase-space distribution function we obtain ~ 30(140) 
/fm 3 as the initial number density, ~ 50(520) GeV/fm 3 as the initial energy density and ~ 
520(930) MeV as the corresponding initial effective temperature at RHIC(LHC). Our minijet 
initial conditions are obtained above j>r = l(2) GeV at RHIC(LHC). As partons below these 
momentum cut-off values might be described by the formation of a classical chromofield, 
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the f(x,p,to) obtained in this paper is the distribution function for the matter part. Initial 
conditions for lower transverse momentum has to be supplied in form of the strength of the 
chromofield or in the form of field energy density in the non-abelian transport equation to 
study production and evolution of the quark-gluon plasma at RHIC and LHC. 
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